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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to generally to data 
matching and, more particularly, to event sequence 
matching. 

2 . Prior Art 

Monitoring a large telecommunication network can 
result in an extensive log of alarms or other events 
of different types that occurred in the system. 
Similar log files may also be produced in mobile 
commerce systems, in web applications, and in mobile 
services. Such logs, or event sequences, generally 
consist of pairs (e,t), where e is an event type and t 
is the occurrence time of the event. 

The sequences of events in a data flow can be, for 
example, sequences of events (alarms) with their 
corresponding occurrence times in a telecommunications 
network. The purpose of finding similar situations in 
these sequences of events, as in many other data analysis 
applications, is to predict events and to understand the 
dynamics of the process producing the sequence. In these 
applications similarity finding can help to customise 
individual services or interfaces by prediction and 
regularities based on previous behaviour. 
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The problem of finding similar situations can be 
described as follows. With a sequence of events 

and a time / and a window width w, find 
another time s, such that the subsequences of S, defined 
as the subsequences of events of S occurring in the half- 
open intervals (t~w,t] and (s-w 9 s] respectively, from here 
on called slices S(t,w) and S(s,w) of S , are similar. The 
slices are sequences of events. The similarity between 
two slices can be defined using an edit distance notion, 
i.e. the distance is defined as the cost of the cheapest 
possible sequence of operations that transforms one slice 
to another. The operations are insertion and deletion of 
an event and moving an event in time and each operation 
has an associated cost. The edit distance can be computed 
using known dynamic programming algorithm. 

Prior art solutions for finding similar situations using 
known dynamic programming algorithms are computationally 
slow in time because of the high computational complexity 
of these algorithms. Furthermore, assigning costs to the 
edit operations is quite problematic as disclosed in 
"Pirjo Moen. Attribute, Event Sequence and Event Type 
Similarity Notions for Data Mining. PhD thesis, 
University of Helsinki, Department of Computer Science, 
Finland, February 2000". In prior art practices, there 
has also been considerable interest in defining intuitive 
and easily computable measures of similarity between 
complex objects and in using abstract similarity notions 
in querying databases as disclosed in: [1] Gautam Das, 
Heikki Mannila and Pirjo Ronkainen, "Similarity of 
attributes by external probes", in Proceedings of the 4 th 



International Conference on Knowledge Discovery and Data 
Mining (KDD) , pages 23-29, 1998; [2] E.-H. Han, G. 

Karypis, V. Kumar and B. Mobasher, "Clustering based on 
association rule hypergraphs" , in Workshop on Research 
Issues on Data Mining and Knowledge Discovery, 1997; [3] 
H. V. Jagadish, A. 0. Mendelzon and T. Milo, "Similarity- 
based queries", in Proceedings of the 14 th Symposium on 
Principles of Database Systems (PODS), pages 36-45, 1995; 
[4] A. J. Knobbe and P. W. Adriaans, "Analyzing binary 
associations", in Proceedings of the 2 nd International 
Conference on Knowledge Discovery and Data Mining (KDD) , 
pages 311-314, 1996; [5] Y. Karov and S. Edelman, 
"Similarity-based word sense disambiguation", in 
Computational Linguistics, 24(l):41-59, 1998; and [7] 
D. A. White and R. Jain, "Algorithms and strategies for 
similarity retrieval", in Technical Report VCL-96-101, 
Visual Computing Laboratory, UC Davis, 1996. 

With ever increasing amounts of information surrounding 
us in our every day life and the numerous applications, 
services, etc., of which the quality relies on data 
processing, faster and more reliable methods for 
information retrieval, and for yielding added value from 
data, are needed to make better or even new applications, 
services, etc. possible. In many fields of applications, 
time series or ordered sets of data are an advantageous 
way of modelling data for many valuable end results. 

Therefore, it is desirable to provide a method and system 
to efficiently analyze large amounts of data. 



SUMMARY OF THE INVENTION 



The foregoing and other problems are overcome, and other 
advantages are realized, in accordance with the presently 
preferred embodiments of these teachings. 

In the present invention, a simple and fast method of 
mapping a sequence of time series data or ordered data 
into points in a k-dimensional Euclidian space is 
introduced. Furthermore it is shown how the mapping can 
be used as pre-processing for finding similar sequences 
in an effective way for retrieving useful information 
from previous sequences of events that resemble a 
specified sequence of events, predicting future events 
based on similarities between previous sequences of 
events and a current sequence of events in a data flow, 
looking for patterns or regularities in data, etc. 

The present invention relates to a method for finding 
similarities in sequences of events in a data flow 
including means for retrieving useful information from 
previous sequences of events that resemble a specified 
query sequence of events. The present invention relates 
also to a method for finding similarities in sequences of 
events in a data flow including means for predicting 
future events in a data flow based on similarities 
between previous sequences of events in the data flow and 
current sequences of events in the data flow. The present 
invention furthermore relates to a method for data 
mining and information retrieval comprising means for 
looking for patterns or regularities in data. 



In accordance with one embodiment of the invention, a 
method for determining similarity between a first event 



set, the first event set comprising a first plurality of 
event types, and a second event set, the second event set 
comprising a second plurality of event types, is 
provided. The method includes steps of mapping an 
observed event set to a multidimensional vector-Q and 
mapping a query event set to a multidimensional query 
vector-q. The next step is a comparison of the vectors 
for a predetermined similarity according to: | |Q-q| | s 
SV, where SV = a predetermined similarity value. 

In accordance with another embodiment of the invention a 
method of finding a query-subset of events within an 
event set is provided. The event set includes a stream of 
ordered events, each ordered event corresponding to an 
event type e from a set of event types E. The method 
includes the steps of dividing the stream of ordered 
events into a plurality of segments and mapping each of 
the plurality of segments to a corresponding dimensional 
segment -vector . The next step maps the query- subset of 
events to a dimensional query-vector; and then compares 
the dimensional query-vector with at least one of the 
dimensional segment-vectors, and as a result of the 
comparison, makes a determination whether the query- 
vector is similar to the compared segment -vector . 

In accordance with one embodiment of the present 
invention, a system for finding a query-set of events 
within a master-set of observed events is provided. The 
system includes a k-dimension random variable generator 
for generating random variables for each event within the 
event set E for each j=l...k dimensions, where k is 
predetermined. The system also includes an observed event 
segmenter for segmenting the master-set of observed 
events to produce a plurality of observed event segments, 
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(dn , Sn-dim, Sim) .« (dhi , Shi-dhm/ s hm ) , where d=segmented 
observed event, and s = a parameter associated with the 
observed event and a weighting function generator. An 
observed event vectorizer for vectorizing each observed 
event segment (d 1 i,Sn...d lrn ,s 1 J... (d hl/ s hl ...d hm , s hm ) according to 
the weighting function provided by the weighting function 
generator and the dimensional random variable 
corresponding to the segmented observed event d and 
dimension k is also provided. In addition, a query event 
vectorizer for vectorizing the query set of events 
(Qei,Qsi...Qe w Qs w , ) according to the weighting function, the 
dimensional random variable corresponding to the query 
event Qe, and dimension k, wherein Qe = a query event, Qs 
= a parameter associated with the query event Qe is also 
provided. In addition, the system includes a comparator 
for comparing the vectorized query-set of events with 
each vectorized observed event segment and generates the 
observed event segment in accordance with the comparison 
results and predetermined similarity factors. 

In accordance with another embodiment of the invention a 
method is disclosed for finding similar situations in 
sequences of events in a data flow, the method including 
the steps of mapping, using a random projection, slices 
of a sequence of events to multidimensional (k) vectors 
and mapping a query sequence of events to a 
multidimensional (k) query vector. The next step searches 
among the event vectors for such multidimensional vectors 
for which the distance from the query vector is small or 
within a predetermined range or value. The last step 
verifies the event segment using event edit techniques. 



The invention is also directed towards a program storage 
device readable by a machine, tangibly embodying a 
program of instructions executable by the machine to 
perform method steps determining similarity between a 
first event set, the first event set comprising a first 
plurality of event types, and a second event set, the 
second event set comprising a second plurality of event 
types . The method includes steps of mapping an observed 
event set to a multidimensional vector-Q and mapping a 
query event set to a multidimensional query vector-q. The 
next step is a comparison of the vectors for a 
predetermined similarity according to: | [Q-q| | s SV, 
where SV = a predetermined similarity value. 



BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing aspects and other features of the present 
invention are explained in the following description, 
taken in connection with the accompanying drawings, 
wherein: 

Figure 1 is a timeline sequence of events A, B , and C 
observed during a twenty second window; 

Figure 2 is multiple timeline sequence of events for 
subsequences, s A ,s B ,s c of sequence s of events shown in 
Figure 1; 

Figures 3A-3C are graphs of two independent event streams 
and their cross-correlation; 

Figures 4A-4D are graphs of an a priori event 
distribution, two independent streams of events following 
the same joint density, and their cross-correlation, 



respectively. The probability of observing 10 or more 
pairs in the analyzed lag intervals is .19% for a uniform 
model, but is 68% considering the a priori event density 
shown in Figure 4A; 

Figures 5A-5B are point graphs showing pairwise 
dependencies for the first ten most frequent events of 
sample telecommunications alarm data; error probability 
threshold is one percent {S= 0.01); Figure 5A shows the 
pairwise dependency with the use of burst information; 
Figure 5B shows the pairwise dependency using only the 
yes/no burst information; dashed lines show weaker 
dependencies; and demanding greater confidence (smaller 
S) renders these events independent ; 

Figures 6A1-6A3 and 6B1-6B3 are graphic examples of 
synthetic data on a small interval; Figures 6A1-6A3 
represent uniform a priori event density, observed 
events, and approximated density (averaging on 30 seconds 
interval) , respectively; Figures 6B1-6B3 represent burst 
type a priori event density, observed events, and 
approximated density (averaging on 30 seconds interval) , 
respectively; note that density scales differ; on graphs 
6B1-6B3 the increase of activity of event #2 after event 
#3 (in about 100 seconds -four last events) and after 
event #4 (at 810,850 and 875, seconds) can be seen; 

Figure 7 is a flow chart showing method steps for 
incorporating event matching features of one embodiment 
of the present invention; 

Figure 8 is system diagram illustrating pattern matching 
features of the present invention; 
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Figures 9 and 9A are system diagrams illustrating pattern 
preconditioning features of the present invention; 

Figure 10 is a scatter diagram of test data showing the 
location of the closest slice to the query slice for 
5 every 10 th event in the test sequence and for W=1000 sec; 

Figure 11 is a collection of 4 scatter diagrams for 
various values of k, showing distances to selected 
windows of example alarm data; 

Figure 12 is a collection of 4 scatter diagrams showing 
10 distances of target windows against the density ratio for 
various values of k for example alarm data; 

Figure 13 is a collection of 4 scatter diagrams showing 
ranks of target windows against the density ratio for 
various values of k for example alarm data; 

15 Figure 14 is a collection of 4 scatter diagrams for 
various values of k, showing distances to selected 
windows of example Entree Chicago data; 

Figure 15 is a collection of 4 scatter diagrams showing 
distances of target windows against the density ratio for 
20 various values of k for example Entree Chicago data; 

Figure 16 is a collection of 4 scatter diagrams showing 
ranks of target windows against the density ratio for 
various values of k for Entree Chicago data; and 

Figure 17 is a table showing distances of fifty closest 
25 windows for k=15, alarm data. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 



The method and system described herein in the present 
invention decreases the amount of needed computations and 
thus computational time by projecting or mapping slices 
to diminish the dimensions of data. The projections being 
such that information of the slice structure is 
preserved, thus making similarity finding possible and 
effective . 

Referring to Fig. 8, there is shown a block diagram of a 
system having an observed event segementer 23 and an 
observed event vectorizer 25 for mapping event segments 
to j -dimensional Euclidean space. The observed event 
vectorizer 2 5 maps the observed event segment to a 
particular Euclidean dimension by selecting a random 
variable for the observed event in a particular dimension 
j from random variable generator 22. In alternate 
embodiments the random variable generator 2 2 may be any 
suitable generator, or it may be a predetermined look-up 
table. Continuing, the observed event vectorizer 25 
weights the random variable for the particular dimension 
with a predetermined function. The predetermined 
weighting function may be any suitable linear or non- 
linear function generated by weighting function generator 
210. The vectorizer 25 continues to weight each random 
variable for each observed event for the particular 
dimension, summing each weighted random variable with the 
previous weighted random variable. Similarly, a query 
event data slice 2 6 is vectorized by query event 
vectorizer 27. The vectorized query event is then 
compared with vectorized observed events by comparator 
28. When comparator 28 detects a similarity within 
predefined limits, the pre-vectorized observed data 
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segment is output. In alternate embodiments the pre- 
vectorized observed data segment may be further tested by 
standard event edit techniques by an optional verifier, 
29. Although the present invention will be described 
with reference to the embodiments and examples shown in 
the drawings, it should be understood that the present 
invention could be embodied in many alternate forms of 
embodiments . 

Referring now to Figure 7, there is shown a method flow 
illustrating steps for performing the event searching 
features shown in Figure 8. Using random projections, 
portions or slices, step 2, of events of S that have 
happened during the last w time units, or a slice of 
width w defined as an sequence of events eeE of 
interest, step 1, are mapped to Euclidean k-dimensional 
space. Mapping, step 4, maps the event slices to k- 
dimensional vectors, step 5 and the query sequence of 
events, step 3, to a query vector, step 7. The method 
then searches, searching, step 6, using a normal linear 
search or more advanced data structures to be discussed 
later, for the vectors, step 8, that are closest to the 
query vector from step 7. The closeness, or similarity, 
meaning that the norm of the difference between found 
vectors, step 8, and the query vector, step 7 is small or 
within a predetermined range or metric. Then, by applying 
edit distance computation for each slice of events 
corresponding to a vector found in step 8, verifying, 
step 9, its similarity to the query sequence from step 3, 
Q, step 10. 
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The random projections are used to compute the k- 
dimensional vector representation r(t i9 w) of each of the 
slices S(t l9 w) and the query sequence Q. To perform the 
random projections, or mappings, step 4, for each event 
type eeE and for each j = l 9 ... 9 k let p(e 9 j)be a normally 
distributed random variable with mean 0 and variance 1. 
Let / be a function from [0 9 w] to [0,l] . Given an event 
sequence 1, S = ((e ]9 t ] ) 9 ... 9 {e n ,t n )) , and a time t, denote by 
the slice of S defined as the subsequence of S 
consisting of those events that happen within the half- 
open interval (t-w 9 t]. Let this sequence consist of the 
events (e u ,t u ) 9 ... 9 (e v9 t v ) = (d l9 s l ) 9 ... 9 (d m9 s m ) , where \<u<v<n and 
m = v-w + l. The random mapping, step 4, of S(t 9 w) to the k- 
dimensional vector r{t 9 w) = (y x (t) 9 ...,y k (t)) can be for example 
defined as 

m 

y } (0 = Z /(' - s < )p( d > > J) ■ EQ . Al 

The function / can be chosen to be linear, for example 
f(x) = x/w, in which case there exists an algorithm to 
compute the vectors r(t.,w) for all the slices S(t l9 w) of 
S . In this embodiment we can write 

y J (.t) = ™\tY J p{d l ,j)-Y J S,p{d i j)) EQ. A2 

>=i ,=i 

For all j = \,...,k and all jo = !,...,« let 
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J) , EQ. A3 

/*=] 

and 

/> 

B j(P)=T J s l P( d ^J) EQ. A4 

Further, let ^(0) = 0 and 5,(0) = 0 for all j. Then A } (p) and 
£,(/?) can be computed in time linear in n , the total 
number of events in the original sequence. We have 

y. (t) = w' ] (t(A j (my A J (h - 1)) - (B J (my B J (h - 1))) , EQ. A5 

showing that given A } and B. and the indices m and h 
corresponding to the interval (t~w,t] f the value ;/ y (/)can be 
computed in constant time. 

Another, non- linear example of the function / could be 

the exponential function /(x)=e~ v . The algorithm in this 
embodiment would be: 

m 

y y (t) = exp(-*)(£ exp(s, )p(d l , j ). EQ . A6 

/=i 

Letting 

C j(P)=H QX V(s i )p(d n j) EQ. A7 

and C,(0)=0, then 

(0 = exp(-0(C y (my A J (h - 1)), EQ . A8 

again showing that after a linear-time pre-processing for 
computing C } (0),...,C J (p) , we can compute each y y (t) in 



13 



linear time. It will be appreciated that in alternate 
embodiments any suitable weighting function may be used. 

Thus, having accumulated, possibly a long sequence S, of 
events eeE and defined a query sequence Q of events of 
a relative time interval of length w , the random 
projection can be used efficiently to find the slices, 
(Figure 7, step 10) of S that are similar to sequence Q. 

In the search phase, step 6, using normal linear search 
or more advanced data structures the vectors 8, 
r(/ ; ,w)cr(/ ; ,w), which are close to the vector g, step 7, 

can be found; closeness meaning that the norm w )~ is 

small or within a predetermined metric. 

In the verifying phase, step 9, after finding the vectors 
from step 8, r(t j9 w), the closeness of the corresponding 
slices S(t j9 w) to the sequence Q can be verified and 

refined, using edit distance computations, yielding the 
final result, step 10. 

Thus, the method presented in the present invention 
avoids the need to do computationally complex full edit 
distance computations in the space of all slices of S by 
doing first the computationally more efficient search in 
the k-dimensional Euclidian space. 

As an example the method described is applied on 
telecommunications alarm data. In one example we consider 
a query slice that can be presented in the form of a 
table as: 
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Table 1 



Event number 


Alarm type 


Relative 
time 


1730 


7260 


0 


1729 


7277 


29 


1728 


1585 


224 


1727 


1940 


821 



Applying the method the present invention describes, the 
following slice is found to be similar with the query 
slice . 



Table 2 



Event number 


Alarm type 


Relative 
time 


29212 


7277 


0 


29211 


7260 


1 


29210 


1585 


205 


29209 


1940 


536 


29208 


7403 


646 


29207 


1903 


825 


29206 


7711 


915 


29205 


7705 


982 



In another example the following query slice shown in 
Table 3 is used to test the method described in the 
present invention : 



Table 3 



Event number 


Alarm type 


Relative 
time 


38000 


1553 


0 


37999 


1553 


0 


37998 


1553 


300 


37997 


1553 


300 


37996 


7002 


313 



15 



37995 


7701 


489 


37994 


7002 


553 


37993 


1553 


600 


37992 


1553 


600 


37991 


1553 


900 


37990 


1553 


900 



Using the method described in the present invention, the 
following slice shown in table 4 is from the sequence of 
events is found to be similar to the above query slice. 



Table 4 



Event number 


Alarm type 


Relative 
time 


37988 


1553 


0 


37987 


1553 


0 


37986 


7002 


72 


37985 


1553 


300 


37984 


1553 


300 


37983 


1553 


600 


37982 


1553 


600 


37981 


7701 


621 


37980 


7002 


692 


37979 


1553 


900 


37978 


1553 


900 



This slice gives the minimum distance at the position 
37988 of the telecommunications alarm data. The distance 
between the slice and the query slice, with k=9, is 
0.107156, computed as the distance between the vector 
representations of the query slice and the found similar 
slice, calculated according to the description of the 
present invention . 
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Known edit distance computations take time 0(nm) for 
slices of lengths n and m , with fairly large constant 
factors. The time and space complexity of the method 
5 described in the present invention is 0(kn) , n being the 
length of the query slice and k being the dimension of 
the used Euclidian space. The larger the length of slices 
is the more is gained in terms of decreasing 
computational complexity when using the method described 
10 in the present invention. In practical embodiments k can 
be chosen to be noticeably smaller than the length of the 
slices without loosing the perseverance requirement of 
the random projections and the method of the present 
invention . 

15 

For the perseverance feature of the method described in 
the present invention, consider the following cases: 

Given two identical slices S(t 9 w) and S(t\w) , the random 
20 vectors r(t 9 w) and r(t\w) are obviously identical and thus 
the distance between the two vectors is zero. 

For two slices that are quite close to each other, 
consider the slices S(t 9 w)=((d {9 t l } 9 (d 29 t 2 )^ t . 9 {d m J m }) and 
25 S{t\w)^{(d 2 J 2 ),... 9 {d m9 t m )) t i.e. the slices are identical 
except that the latter does not have the element . 
For simplicity, assume that f(t-t l ) = l for all t i . Then the 
squared distance between vectors r(/,w) and r(t\w) is 

k 

d 2 =^r(d n if . Since the variables p{d x j) are independent 

30 identically distributed, the expected value of the 
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squared distance is E(d 2 )= kE(p(d ]9 \f) = k . The variance is 
E(d A )-E{d 2 ) 2 =k{3-\ 2 ) = 2k. 



In contrast, consider the squared distance from the 
origin to a vector obtained by mapping a random m-event 
slice S(t,w) = ({d x ,/,),..., (d m ,t m )) . Every element of the vector 

r(t,w) is a sum of random variables, r.(t 9 w) = ~t J )p(d J ,i) . 

Assuming they are independent, i.e., that no event is 
repeated in the slice, the expected value of r t (t 9 w) 2 is 

E (n(t> w ) 2 )=i E [p( d so the expected squared 

distance from the origin to the vector is km . For the 
variance, we get 

4;M 4 )-4,M 2 ) 2 =t[E{p(d J 4 4 -E( P (d j , i y) 2 }= m [3-l^2rn , 

7=1 

so the variance of the squared distance is 2km . Thus the 
expected difference between the two similar slice vector 
representations relative to the expected length of the 

vectors is of order (m)~2 , i.e. for practical values of m 
the expected distance between vector representations of 
the similar slices is small. 

Thus the editing operations of inserting and deleting one 
event have small expected effect on the distance, 
compared to arbitrary vectors in the k-dimensional space. 
In the previous analysis, we assumed that all the events 
have equal weight. In practice, the effects of these 
editing operations will be even smaller for events 
towards the low-weight end of the slice. Also, assuming a 
continuous function /, the editing operation of moving an 
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element in time has an effect proportional to the length 
of the time change, just as in the notion of editing 
distance. Therefore, the approximated distance is small 
when a slice is subjected to a small number of editing 
operations . 

The inverse relationship can also be shown. That is, if 
two slices (i.e., segments) are far from each other, then 
the corresponding random vectors are far from each other 
with high probability. 

Examples of the present invention follow. First, recall 
that the goal is to locate previously occurring 
situations that resemble the current one. To test how 
well this goal is met, the following experiments for 
alarm data are shown. The alarm data consists of 46662 
alarms over a period of a month. The time values seen on 
the x-axis of several figures are in seconds, and they 
range from about half a million to three million. The 
figures show only a subset of the full range for reasons 
of legibility. There are 180 types of alarms occurring in 
the data. The Entree Chicago data comes from the log of 
an application which provides restaurant recommendations 
on the Web. Here, we consider only the final 
recommendation made by the system for a single user. The 
data is converted from the date/ time information to 
seconds since the beginning of 1970 (a popular way of 
expressing time in some operating systems , so the time 
values vary from about 842 million to 924 million, and 
again, the figures show a smaller range. The sequence 
contains 50672 events of 619 types. 
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First, to get a qualitative sense for the alarm test 
data, all closest slices for every 10 th window with w = 
100 0 using the random mapping approximation was 
calculated. Figure 10 shows the location for the closest 
5 match for each slice. Some temporal locality can be 
observed. Next, an artificial query window of width 1000 
consisting of 37 events of 20 types was created. None of 
the event types occurred in the original sequence and 50 
copies of this query window were pasted onto the alarm 

10 data. The pasting was made transparently, i.e., the 
original events were left into the sequence. Thus the 
modified sequence contains some clear copies of the query 
window, but most of the copies contain also other events. 
The query window was constructed by copying a part of the 

15 sequence and renaming the events. Next, exercising the 
mapping features of the present invention, with several 
values of the parameter k, and performing queries with 
the following rule: list all windows in the order of 
approximate distance to the query window, except when a 

20 window S is listed, disregard all subsequent windows that 
overlap S. For example, if (9000, 10000] is listed as one 
of the resulting windows, none of the windows {9000+i, 
10000 + i] for i = 999,... , 999 are listed. 

25 Figure 11 shows the distances of some of the resulting 
windows for k = 5, 10, 15, 30. The x-axis gives the 
location in the sequence, and the y-axis gives the 
approximated distance, i.e., the Euclidean distance in k- 
space. The plots have been truncated at a point slightly 

30 above the median distance for scaling purposes, since 
some windows have extremely high distances. For 
legibility, the horizontal axis shows only a art of the 
full time scale. The vertical lines indicate the 



positions of the inserted windows. The best 50 non- 
overlapping windows in the case k = 15 are listed in the 
table shown in Figure 17. The first column gives the 
computed distance to the query window; in the first three 
5 cases, there were no extra events intermingled with the 
target, so the distance was zero. The second column 
indicates the position (i.e., time in seconds of the 
window, and the third one gives the position of the 
closest target window, if there is one within the window 
10 width 1000. The fourth column is simply the difference of 
the second and third ones. Of the 50 target windows, 22 
M= are found exactly, and 20 more overlap one of the best 50 

S windows. Note how all the distances below 13 correspond 

W to windows close to the targets, and distances up to 

m 15 about 5 are all exact hits. This is illustrated in Figure 
12, where the approximated distances of the inserted 
s target windows are plotted against the "density ratio" 

|l'|/|l| , where | I | is the number of inserted events 
yj and | I' | is the actual number of events within the 

!Z 20 1000-width target window, counting both inserted events 
M= and those that were art of the original sequence. 

For each target window W, either W appears or one or two 
windows overlapping W appear in the sorted list where 

25 overlaps have been omitted. We call the "rank" of W in 
this listing the rank of the first listed window W ' r\W'+# . 
Ties are resolved arbitrarily to give each window a 
unique rank. Figure 13 shows the ranks of target windows 
against the density ratio defined earlier. Ranks greater 

30 than 100 are shown as 100. Note the effect of kz with a 
larger number of dimensions, the targets stand out better 
in Figure 11, but the number of targets within the best 
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windows doesn't increase very much. With this data, it 
seems that values of k around 5 orlO are sufficient. 

In the case of the Entree Chicago data, we again inserted 
50 windows in the sequence. Since the time scale is 
different, we used 10000 seconds as the window width. In 
contrast to the previous experiments, the inserts contain 
mostly event types that also appear in the original 
sequence: of the 24 events in the insert, 10 are 
artificial (ones that don't appear in the unmodified 
sequence. This makes the targets a bit harder to 
recognize, but they should still be quite far away from 
other windows. As in the previous case, the insert was 
obtained by taking a art of the sequence and renaming 
some events. Figure 14 shows the distances of windows 
selected as in the case of alarm data. The values of k 
shown are 10, 15, 20 and 30. The lots indicate that 
several of the targets can be found in this case as well. 
With this data and this insert, the density ratio of the 
targets was higher than in the previous case. This is 
nicely shown in Figures 15 and 16, plotted for the cases 
k =10and7c = 30 . 

Referring to Figure 9, an alternate embodiment may 
decompose the observed and query events into independent 
and dependent events before vectorizing the event 
segments. The decomposability of an event sequence may be 
stated as whether there are two (or more) unrelated sets 
of event types such that the corresponding event 
subsequences are independent. In this embodiment, an 
observed event data pre-conditioner , 42, preconditions 
the data into independent event segments 43 and dependent 
event segments 44. Similarly, the query event segment is 
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separated by query event data pre -conditioner to segment 
the query event slice into independent query events and 
dependent query events. Vectorizers, 25A,25B and 27A,27B, 
operate as described above, as do comparators 2 8A and 
28B. 

Referring to Figure 9A, an alternate embodiment may 
precondition the observed and query event sets 9A41 and 
9A26 according to a statistical analysis of the query 
event set 9A26 by query event data pre-conditioner 9A46. 
The output of query event data pre-conditioner 9A4 6 is 
then used as an input to observed event data pre- 
conditioner 9A42. It will be appreciated that in 
alternate embodiments, observed event data pre- 
conditioner 9A46 may statistically analyze the observed 
data set 9A41 and use the results in conjunction with, or 
independent of, the statistical analysis results provided 
by query event data pre-conditioner 9A46. Vectorizers, 
9A25 and 9A27 operate as described above, as does 
comparators 9A2 8. 

Referring again to Figure 9, the event data pre- 
conditioned for separating events into independent and 
dependent sets are described as follows. 
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We consider the input as a sequence of events, where each event luw an 
associated time of occurrence. Given a set E = {ei,... ,c*} of event types, 
an event is a pair (A, i), where A G is an event type and I 6 N is the 
occurrence time of the event. Note, that we often use the term event referring 
to the event type; the exact meaning should be clear from the context. 
An event sequence s on E is an ordered sequence of events, 

s = ((A il t l ),(A 2 ,t 2 ),... ) Un,t n )) (i) 

such that At € E for all i = 1, . . . ,n, and U £ [T s ,T e ], ti < for all 
i = 1, . . . , n — 1, where T s , T e are integers denoting the starting and ending 
time of the observation. Note that we can have t x — t 1+ i, i.e., several events 
can occur at the same time. However, we assume thafr for any A 6 B at 
most one event of type A occurs at any given time. 

Given an event sequence $ over a set of event types B, and a subset 
Ei C E, the projection s[Ei] of 5 to E\ is the event sequence consisting 
of those events (e, t) from 5 such that e € E\. A sub-sequence of event e i5 
denoted by s Bi , is a subsequence of s consisting only of the events of type e x 
from s, i.e., s e , is a projection of 5 onto E\ — {ei}. 

Alternatively, we can view s as a function from the observed period, 
[T 5 ,T e ], into {0,1} |£| , and {s e .} eieE as functions from [T s ,T e ] into {0,1}, 
such that $ = s ei x ... x s efc . In such formulation, s(t) denotes the events 
that happened in the time unit t. 

EXAMPLE : Figure 1 presents the event sequence of three event types E = 
{A, B, C} observed for 20 seconds, that is T s = 1, T e = 20 and 

a = <(B, 1), (C, 2), (A, 3), (A, 5), (A, 8), . . . , (B y 20), (£7, 20)). 

Note that a number of events of different types can occur in the same second. 
The subsequences of sequence s are shown on Figure 2 and they are 

s A = ((A 3), (A, 5), (A 8),..., (A, 18)) 

* B = {(^, 1), (B, 9), (B, 13), (B, 18), (B, 20)) 

3c=((C,2),(C J ll),(C,14),(C,20)) 

It can be seen that event C always follows event B with one or two sec- 
onds lag. The C event that follows (B, 20) was not observed due to finite 
observation time. 

Treating 5 as a function from [1,20] into {0,1} 3 we have 

s = (010, 001, 100, 000, 100, . . . , 000, 011) 

and sa, sb and sc are just a binary vectors of length 20: 

s A = 00101001000111010100 
s B = 10000000100010000101 
s c = 01000000001001000001. □ 



Decomposition of event sequences 

In order to discuss the independence properties we are interested in, we have 
to provide a way of probabilistic modeling of event sequences. 

Given a set E of event types, the set of all event sequences over E can 
be viewed as the set T E of all the functions Z : [T S ,T € ] -> {0, 1} |K 1 That is, 
given a time £, the value Z(t) indicates which events occur at that time. 

A probabilistic model for event sequences is, in utmost generality, just 
a probability distribution \i E on Te> For example, y. E may depend only 
on the total number of the observed events and give a higher probability to 
the sequences that contain N events, i.e. fi E (Z) = a • e "^ N ~ Nz ^^ 2 where 
N z denotes the total number of events in Z, N z = Y^iLr, ll z WIUt and 
a, b are some appropriate constants. Note that in this example all event 
subsequences are dependent. 

Next we define what it means that a distribution of event sequences is an 
independent composition of two distributions. We use the analogous concept 
from the distribution of discrete random variables: Let {Xi,... y X p } be a 
discrete variables and denote by P(X L = x u . . . ,X p = x p ) the probability 
of observing the value combinations (i lr .,,i p ). Now P is an independent 



composition of distributions over variables {X u . . . , Xj) and {Xj h, • ■ • i 
if for all combinations . . . ,z p ) we have 

P(Xx = x l5 . . . ,X p = x p ) = 

(2) 

where Pi and P 2 are the marginal distributions defined by 

Pi(Xi = x u . . . ,X 3 = xj) = 

Y.(x : + ly ...,x p ) P ( X l -X U ...,Xj =Xj,Xj+i =X3 + u...,X p = x p ) 
P L {Xj+i s=Xj,... 9 X p = x p ) = 

D{xi ,...,*,) ^(-^i =xi».-.,-X' (7 - =^,^+1 s= £;+!,..., ^ =a:p)* 

* (3) 
The above definition is, of course, easily extended for the decomposition of 
{X\ , . . . , X p } into more than two subsets. 

Now, let Ei be a subset of E. The distribution he defines naturally the 
marginal distribution liEi on Te x 

MEi(*l) = ^ ^ 

s€^£r,s[^i]=5i 

We can now provide a decomposition definition: 

Definition [Event set decomposition]: The set of event types E 
decomposes into pairwise disjoint sets E\, . . . ,E m with E = U t =i 
and Vi ^ j, P z 0 P 7 = 0 t/ /or a/J 5 € J^; 

m 

T/iat is, the probability of observing a sequence s is the product 
of the marginal probabilities of observing the projected sequences 

s[E,}. 

If E decomposes into E u . . . , E m , we also say that ue decomposes into 
PBi » • ■ • > PEm aiK * ^at & consists of independent components J5?i , — , Em- 
As a special case, if E consists of two event types A and P, it decomposes 
into A and B provided 

I.e., the occurrence probability of a sequence of A's and P's is the product 
of the probability of seeing the A's and probability of seeing the P's. 
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Finding independent components from observed sequences 

Our goal is to start from observed sequence s over a set of event types E 
and to find sets Ei,...,E m such that the probability distribution fi E on 
Te is decomposed into the marginal distributions /U£ t , . . . ,A*E m - There are 
two obstacles to this approach: First, we only observe a sequence, not /Z£. 
Second, the set of alternatives for E\, . . . , E m is exponential in size. 

The first obstacle is considered in Section 3.1 where we show that certain 
quite natural conditions can be used to obtain information about fj,£ from 
a single (long) sequence over E. 

We next describe how to cope with the second obstacle. We overcome 
this problem by restricting our attention to pairwise interaction between 
event types. That is, given fjt E , two event types A and B are independent, 
if for all s £ Te we have 

t*lA,BM{A,B}]) = pa(sa)-I*b{sb). (?) 

We show in the next section how we can effectively test this condition. 

Given information about the pairwise dependencies between event types, 
we search for independent sets of event types. Let G = (E, H) be a graph 
of E such that there is an edge between event types A and B if and only 
if A and B are dependent. Then our task is simply to find the connected 
components of G, which can be done in 0(\E\ 2 ) by any standard algorithm 

Using the above procedure we separate E into the maximal number 
subsets E U ...,E U such that VI < i ^ j < /,Ve' 6 £»,Ve" 6 E 3 : e',e" are 
independent. Note, that pairwise independence generally does not imply the 
mutual independence^ In our case it means that E u ... : E t 

is not necessarily a decomposition of E. We use, however, E\,.>.,Ei as a 
practical alternative to a true decomposition of E. In the remainder of this 
paper we will concentrate on detecting pairwise dependencies among the 
events. 

Detection of Pairwise Dependencies 

The definition of decomposability given in the previous section is based on 
the use of the distribution fi E on the set of all event sequences. This makes 
it impossible to study decomposability of a single sequence. If we have a 
large set of observed sequences, we can form an approximation of p^. Given 
a sufficiently long single sequence we can also obtain information about p£. 
In the following subsection we describe the conditions under which this is 
the case. 
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We expand our definitions a bit. Instead of considering event sequences over 
the finite interval [T 5 ,T e ] of time, we (for a short while) consider infinitely 
long sequences. Such sequence 5 is a function Z {0, 1} |£;I , and s(t) gives 
the events that happened at time t. 

We assume that the event sequence is generated by some underlying 
stochastic process {Z t } te z, where Z t is a random variable that takes values 
from {0, l} j£| . In this formulation T E is a set of functions from Z into 
{0, 1} |£;| , T E = {Z{t)\Z(t) : Z -> {0, and pt E is a probability measure 
on Thus, the observed event sequence 5 is some specific realization 
f(t) e T E restricted to the interval [T s1 T e ]. 

First two assumptions that we introduce will permit to draw general 
conclusions from the single log, while the third assumption will allow us to 
restrict our attention to the local properties of the event generation process. 

Assumption 1 (Stationary Process) The observed process is a station- 
ary process, i.e., it is shift-independent: 

Pe(S) = me(5 +t ), Vt€Z,V5C^b (8) 

where S +T = {f+r(t)\3f <= S,s.t.Vt e Z : f +r {t) = f{t + r)}. 

The assumption of stationary process means that process does not change 
over time. While this assumption by itself is somewhat unrealistic, in prac- 
tice it can be easily justified by windowing, i.e., considering only a fixed 
sufficiently large time period. The question of stationary testing for a spe- 
cific stochastic process is of great interest by itself, but it is beyond the scope 
of this paper. 

Assumption 2 (Ergodic Process) The observed process is an ergodic pro- 
cess, i.e., statistics that do not depend on the time are constant. That is, 
such statistics do not depend on the realization of the process. 

This is a very important assumption that means that any realization of the 
process is a representative of all possible runs. In particular it means that 
we can average by time instead of averaging different runs of the process. 
. Let X{f,u) denote the time average 

fT 

X(f> u ) = J™ (l/T) / f(u + t)dt (9) 

i — ►oo J —T 

where / 6 Te is a particular realization (event-log). This random variable 
is time invariant. If the process is ergodic, then X is the same for all /, i.e., 
X{f,u) = X, and for a stationary process we have 

X = E[X(f, u)} = lim (l/T) f T E[f(u + t))dt = / (10) 



where / = f(t) = E[f(t)], so the expected value in every point, /, is equal 
to the time average X. 

Note that not every stationary process is ergodic. For example, a process 
that is constant in time is stationary, but it is not ergodic, since different 
realizations may bring different constant values. 4 

The assumption of ergodicity is very intuitive in many natural systems, 
e.g., in telecommunications alarms monitoring. In such systems, we feel that 
logs from different periods are independent and are a good representative of 
the overall behavior of the system. This observation is also the basis for the 
next assumption. 

Assumption 3 (Quasi-Markovian Process) The observed process isquasi- 
Markovian in the sense that local distributions are completely determined by 
the process values in some finite neighborhood, i.e. 

P (Z t e D\z f ,t f ^t)= P (z t e D\z t ,,t' ^t,\t- 1'\ < K) (11) 

where D C {0,1} |£:| and K is some predefined positive constant, which is 
called maximal lag. 

We call this assumption Quasi-Markovian in order to distinguish it from 
the classical definition of Markovian process where K = 1. We specify 
that local probabilities depend not only on the past, but also on the future 
to account for cases with lagged alarms and alarms that originate from 
unobserved joint source but have variable delay times. 

Note that Markovian property does not say that random variables that 
are too far apart (i.e., lagged by more than K second) are independent. 
It simply says that the information that governs the distribution of some 
particular random variable is contained in its neighborhood, i.e., in order for 
one variable to have an influence on another over the maximum lag period 
this variable should "p&ss^ the influence information in time steps smaller 
than K seconds. 

First order dependencies 

The straightforward way to detect pairwise dependencies among the events 
is by direct test of the pairwise independence condition. However, such 
approach is infeasible even for the simplest cases: Consider that two events 
are generated by stationary, ergodic and quasi- Markovian process with AC = 
30 seconds. In this case, we would like to approximate probabilities of the 
event distribution on some arbitrary 30 seconds interval (the start- time of 
the interval is unimportant since the process is stationary). This task will 
require approximation of probability of 2 30 -2 30 « 10 12 joint event sequences. 
Supposing that the average of 100 observations of each sequence are needed 



to approximate its true frequency one should observe the event generation 
process for about 10 14 seconds, which is approximately 31 million years. 

The example given above demonstrates that there is no feasible way to 
detect all possible event dependencies for arbitrary event generation process. 
For many inter-event dependencies, however, there is no need to compute 
the full probabilities of event distribution functions on interval £, since the 
dependencies among the events are much more straightforward and are de- 
tectable by simpler techniques. For example, one event may always follow 
another event after a few seconds (see example on Figures 1,2). Such de- 
pendency, called episode, is easily detectable. 

This work deals with detection of event dependencies of first order. Such 
event dependencies can be described by specifying the expected density of 
events of one type in the neighborhood of events of second type. These 
neighborhood densities can usually be approximated with sufficient preci- 
sion given the typical number of events (hundreds) in the data streams that 
we have encountered. Note also, that in the many applications the event 
streams are very sparse so it is reasonable to calculate densities in the neigh- 
borhood of events and not in the neighborhood of 'holes 1 (periods with no 
events occurring). Otherwise, the meaning of event and not-event may be 
switched. 

Cross-correlation analysis 

Consider two events ei and e2. We observe a joint stochastic process that 
consists of two (possibly dependent) processes: one is generating events of 
type ei and second is generating events of type e2- Consequently we have 
two streams of events si, s 2 of first and second event respectively. We can 
view si and 52 as a functions from the observed time period [1;T] (where 
T is the length of observation) into event frequencies, {0,1}. An example 
of such process is given on Figure 3. 

Supposing the quasi-Markovian property of the event generation process, 
the first order dependency should expose itself in the 2K + 1 neighborhood 
of each event. We define the cross correlation with maximum lag K, and 
with no normalization: 



Note that the cross correlation vector C12 is the reverse of C21. By dividing 
C12 by the observed frequencies of ei and e2 we get the estimate of the 
neighborhood densities of e2 in the neighborhood of ei and of e\ in the 
neighborhood of e%. 

Ideally, if two events are unrelated and the length of observation (T) is 
sufficiently large, the average density in the event neighborhood should be 




m > 0 
m < 0 




(12) 



the same as everywhere on the observed period. It is the same as to require 
that lagged cross-covariance is everywhere zero, i.e. 

cov l2 {m) = ci 2 {m)/(T - m) - pip 2 = 0, Vm € [-/C, /C], (13) 

where pi, P2 are the a priori event probabilities, that does not depend on 
the time of observation since the process is supposed to be stationary. These 
probabilities can be estimated by averaging the observed frequencies of e\ 
and e 2 over the length of the process (this is the direct usage of ergodicity 
assumption), i.e. let rji, n 2 denote the observed number of events ei and e 2 
respectively, thus 

Pi^m/T (14) 

In practice, the neighborhood densities are deviating from the expected 
values even if events are independent; this is due to the random nature of the 
event generation process and due to finite number of observations. Thus, we 
should introduce some model that will account for these effects and give us 
a threshold values, that will allow detection of the event dependencies that 
are beyond random phenomenon. 

Modeling the independent event generation processes 

Consider two independent, stationary stochastic processes that are generat- 
ing events of types e\ and e 2 . We assume that the each event generation 
process is not auto correlated, i.e., in each process the probability of event(s) 
occurring at any given time is independent on the nearby events. Such as- 
sumption maybe justified in the case of sparse, quasi-Markovian processes 
where the average distance between the events of the same type is large 
comparing to the maximum lag distance. 

We are interested in computing the probability of encountering c\2{m) = 
k for some particular m over the observed stream of length T. Since the 
event generation processes are assumed to be independent and stationary 
the above question is equivalent to calculating the probability of observing 
c l2 (0) = A:. We are also interested in computing the probability that c\ 2 {m) 
will not exceed some predefined values on the range m € [— K,K]. We 
formulate the following problem: 

Problem 1 (Pairs Distribution) Assume we observe 2T independent bi- 
nary random variables s\(t), s 2 {t) for t = 1,. . . ,T, with P($i(t) = 1) = Pi 
and P{s 2 {t) = 1) = P2 for all t — 1, . . . ,f. Let c\ 2 be defined by Equation 12. 
The questions are: 

• What is the distribution of Ci 2 (0)? 

• What is the exact from of P(Vm € [A,B],a < c i2 (m) < b) for some 
given {A,B,a,b}? 
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We answer the first question exactly and give the approximation scheme 
for the second. 

Under the assumptions of Problem 1 the generation of pair of events 
e\ and e 2 is independent on its neighborhood and the probability of events 
et, e2 occurring together is p^ = piP2> where pi, p2 are a priori event 
probabilities. Thus the probability of observing exactly k pairs of events ei, 
e 2 is described by binomial distribution: 

P(c, 2 (0) = k) = (J)p? 2 (l -Pi2) T - k (15) 

To assess the probability of observing a random phenomenon we would 
like to estimate the probability of obtaining more or equally extreme number 
of observations as k, i.e. 

11 = P(random phenomenon) = P(|ci 2 (0) - T-pn] > \k - T • P12I) (16) 

Direct calculation of 1Z may be a hard computational task, which is un- 
necessary since we can use one of the known approximations to binomial 
distribution, namely to approximate binomial distribution by Normal or 
Poisson distributions. Since the typical values of pn encountered in prac- 
tice are very small (for the two most frequent events in telecommunications 
alarms datapi2 = 1.6 -10~ 6 ) the Poisson approximation is more appropriate: 

p( Cl2 (o) = k) « u = r- Pl2 . (17) 

Thus the risk of accepting a random phenomenon with lag m as a true 
correlation is 

7l(c 12 (m) = fc)«^f^p«l_^^ (18) 

The minimal extreme values are not considered by the above formula, since 
for the typical data we worked with the probability of observing zero pairs 
is quite large. 1 It is worth to mention, that we have only observed T — m 
trials of the events lagged by m seconds, but T is usually much bigger than 
m so the error of assuming T trials is insignificant . 

We approximate the probability of 012(771) to achieve the certain values on 
the range of by assuming that c i2 (m) and ci 2 (m') are independent 

for m ^ ml. We have: 

P(Vm €[-£,£}, a <ci 2 ( m) < b) ^ [P(a < c I2 (0) < b)] 2fC+l . (19) 

l The usage of particular approximation (Poisson or Normal), as well as, the usage 
symmetric or asymmetric risk calculations are dictated by particular application. 
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Lot. S denote the tolerable error probability of accepting two events ;is de- 
pendent, so threshold on the risk of each particular lagged correlation should 
be set as: 

EXAMPLE: Consider two independent event streams shown on Figure 3. We 
would like to know what is the random probability of the observed event: 
'The maximum number of 7 pairs was encountered while analyzing lagged 
pairs, for lags in [-20, 20]\ Under the model described above the probability 
of observing 7 or more pairs for one particular lag is (assuming lag is much 
smaller than observation time): 

P{#pairs > 7|lag = - 17) « V ~- = 0.0335. (21) 

Thus, assuming the pair generation trials were independent, we have: 

P{#pairs > 7j - 20 < lag < 20) = 1 - (1 - 0.0335) 41 = 0.7528. (22) 

So the probability of observing 7 (or more) pairs in the analyzed lag in- 
terval [-20,20] is 75% for data in Figure 3, thus these streams can not be 
considered dependent. □ 



In-burst event independence analysis 

In many event sequences events tend to appear in bursts. Burst is a sudden 
increase of activity of the event generating process. For example in the 
telecommunication alarms data that we have analyzed, the mean inter-event 
distance is 68 seconds, the median is only 20 and maximal time between 
subsequent events is 3600 seconds(l). This data indicates that alarms in the 
telecommunication network data tend to appear in bursts with long intervals 
between them. 

In burst-type event sequences, most of the events are dependent just 
because they are almost always grouped together. However, we may still 
want to perform in-burst analysis of event independence. Such analysis 
can be seen as deciding on the event independence given that events are 
grouped in bursts. Note that this model describes the situation when bursts 
are ignited by some external event and knowledge of these external events 
may rend many of the in-burst events independent. 

To assess the risk of assuming the random phenomenon as true in-burst 
event dependence we would like to address the following problem, which is 
the generalization of Problem 1: 

Problem 2 (Pairs Distribution, Non-Uniform density) Letrji, n 2 be 
a positive constants and let V(t) : {1, . . . , T} -+ [0, max( * 1>7?2j ] be a function 

with integral one t i.e., X^ t 2>(0 = 1- We observe 2T independent binary 
random variables si(t), $ 2 {t) for t = 1,... } T, such that P{si{t) = 1) = 

rnV{t) and P(s 2 (t) = 1) = ifeD(t) for a//Ul T. Let c L2 be defined by 

Equation 12. The questions are: 

• What is the distribution of c l2 (m)? 

• How to estimate P(Vm € [A, J3], a < c i2 (m) < b) for some {A, B, a, 6} ? 

This problem is illustrated on Figure 4, where the two independent event 
streams are generated according to some a priori density V. Formally speak- 
ing, we assume two-stage event-generation process. First a subset Q oiT ei ,e 2 
is chosen such that expected event density is equal to V y i.e. Ejz${f Ci ) = n{D 
(for i = 1,2), and only then specific / 6 Q is chosen. In this way the ob- 
served event streams are 'forced' to comply with some a priori event density 
V. We would like to find if c u e 2 arc independent given V. 



. Simple in-burst event dependence analysis scheme 

The problem with the above approach lies in the estimation of a priori 
event density, which is too biased to the current realization of the random 

densUv '^Z^ 1° ° VerC ° me this difficult * a » d » ltrod "<* * <°'>'*'' 

e.g., mixture of Gaussian* i 

a priori distribution Te L P fs^pUon is a 'binary' form ol 

information about bursts *' V ^ W 

A r, u j nd m - b urst event density is uniform. 

stream ts vervt ^ * *** S6Cti ° n is based on the fact that evcllts 

niumt W«k AT ParSe V and there *** 1°«S intervals between subse- 

Ma t ny ° f the interv ^ are greater than the maximal lag time, 
dent evTn K ^ "» be Safe ^ *P«*«i a ««»*« of indepen- 
d!at f L f qU6nCeS COTres P°nd to bursts and inter-burst intervals 
that are free from any events. 

n.riISnT dard ana ' ysis ° f • " » Performed on the burst 

Lis ofln^'r ' theseriesof ^nts that are separated by no-event inter- 
° . a V east - Such analysis allows detecting first- order indepen- 
Tl JrZ g T bUrStS (bUt a " mni »B nothin S ^ut burst magnitude), 
onlv " ve mS 6Vent Polities Pl and V , from bursts areas 

of he Znd f ' ^ 63 Pi ^ » ^ Uation 14 ) rendering more 
of the found dependencies random comparative to the 'uniform' analysis 

is a^lfKT m , ^ SUCh SimpHfied in - burst event independence analvsis 
is outlmed below (step one is the same for all event pairs): 

1. Remove all no-event intervals that are greater that JC seconds. 

2. Calculate the lower limit on the number of observed lagged pairs for 
a given J (Equation 18 and Equation 20) using the a priori event 
probab.h ties est.mated from remaining data (Equation 14 with T new = 
T - time(no - event - intervals)). 
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Calculate the observed number of pairs for each lag in [-£,£] and 
announce the events as dependent if the observed number of pairs for 
some lag exceeds lower limit calculated on Step 2. 

An example of in-burst independence analysis 

formed™ C °T id 7 7° ^ StfeamS ° n Fi * ure 4 - Performing the uni- 
diSbtl of C f C (SeCti ° n 3 - 4) With0ut takin S ^ ^ount the 

0 r h C o T* ^ °T Pr0babilHy ° f °- 19% t0 obse ™ cor«Iation 
onlv wi h abo t ,0 0VCr , ha r nd ' rem ° VinS the n °- event ««. ■*» —king 
U^o^ b ^> - - 1* Polity thai 

This analysis shows that two events are clearly dependent in general 
without ons.der.ng an information about prior distribution. The in-burs 
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Note, that this result is achieved under the very simple model, without 
even taking into account the actual form of event density. A natural ex- 
tension may be to threshold the estimated a-priori density function at some 
over label (and not at zero, like in the above approach). This method will 
allow gradual adjustment of the event dependencies, from the events inde- 
pendent regardless to bursts to the events that are dependent even given 
the information that they occur together in some dense bursts (applying 
threshold on the estimated event density function at some high level). 

Experimental Results 

The empirical analysis of the proposed dependency detection methods was 
performed on the telecommunications alarm log data and on two synthetic 
datasets that were especially created to test the dependency detection algo- 
rithm. The data was analyzed using four dependency detection dependency 
detection methods, as summarized below: 

• Event dependency analysis under uniform event generation model, 

• Event dependency analysis using only the yes/no burst information, 

All algorithms were applied with maximum lag K = 300, and error proba- 
bility threshold 8 = 0.01. Parameters of the episode detection algorithm are 
described in the corresponding sections. 

Telecommunication alarms data 

The telecommunications alarm log consists of 46662 alarms in telecommuni- 
cation network logged over the period of one month (2626146 seconds). The 
alarms are of 180 different types and 27 alarms are occurring with relative 
frequency of more than one percent. The original data contains a number 
multiple events of the same time occurring in the same second. We suppress 
these multiple entries to allow only one event of each type in any particular 
second. This operation leaves 38340 events in the log, which correspond to 
82 percent of the original data. The mean inter-event distance is about 68 
seconds, while the median is only 20 indicating that events tend to group 
together in bursts. 

We restrict our analysis to the ten most frequent alarms that are respon- 
sible for more than 51 percent of the alarm data. These alarms, occurring 
more than a thousand times each, are frequent enough to allow various 
meaningful statistical approximations. On the other hand, such restriction 
enables to follow algorithm behavior in detail and not to be overwhelmed 
by large amount of data and inter event relationships. 

We perform analysis of pairwise event dependencies in four stages going 
from the most extensive set of dependencies detectable under no burst as- 
sumptions to the strongest and exclusive dependencies called episodes. The 
results are illustrated in Figure 5. 

We first analyze the telecommunications alarm data using the 'uniform' 
event generation model, i.e., without respect to the burst information. We 
demand random probability of the observed lagged correlation to be less than 
one percent and we are analyzing less than 5 minutes lags , i.e., 5 = 0.01, 
XT = 300 and K th = 1 - (1 - O.OljVMi % L6 7 - 10~ 5 . Pairwise event 
dependencies that were detected in this analysis are shown on Figure 5a. 

Removing the no-events intervals that are longer than JC = 300 seconds 
and applying the dependencies detection technique with the same param- 
eters we get fewer dependencies, as illustrated on Figure 5b. The depen- 
dencies that are dropped are dependencies in pairs (2,6), (6,8) and, most 
noticeable, (1,2) and (4, 10). Note that every in-burst dependency is also a 
dependency in the general sense. 

Note that the set of inter-event dependencies consistently decreases as 
we restrict the definition of dependent events. Figures 5(a,b). 



Experiments with synthetic data 

Wo conduct additional experiments with two sets of synthetic data. The syn- 
thetic data streams contain events of five different types occurring over the 
period of 600000 seconds. The events frequencies are about 5000,3000,1000,1000 
and 100 in the observed streams. Two rules were applied during the gener- 
ation of event streams: 

• Event #2 follows event #3 with 80% probability in time frame [05, 105] 
seconds. 

• Event #2 follows event #4 with 90% probability in time frame [8, 12] 
seconds. 

All other events were generated independently. 

The first data stream was generated with events appear uniformly over 
the whole time period, while the second data stream was generated accord- 
ing to a priori distribution, which consisted of uniform distribution plus 
200 randomly located Gaussians with variances varying from 300 to 3600 
seconds. These datasets are illustrated on Figure 6. 

The same pairwise dependency analysis techniques were applied on these 
synthetic datasets to test the performance and stability of the proposed 
dependency analysis methods. To support the consistency of the results 
the techniques were applied with exactly the same parameters as for the 
telecommunication alarms data, namely 8 = 0.01, JC = 300 and density 
approximation by 30 seconds average of the observed event stream. The 
experimental results are shown in Table 3* In the first dataset, the only 
dependencies in pairs (2,3) and (2,4) are correctly detected by all methods, 
and all other pairs were announced independent. 

In the second dataset, ignoring the burst information renders all of events 
to be dependent, and it is expected, since all the events are inherently depen- 
dent because they follow the same a priori distribution. The simple in-burst 
independence analysis, which takes only burst existence into account, also 
announces almost all (except (5,1), (5,2) and (5,3)} of the events to be 
dependent. Possible source to that behavior can be even higher (comparing 
to telecommunications data) in-burst event density, since the median inter- 
event distance is only 2 seconds, comparing with 20 for telecommunications 
data, while the mean inter-event distances are about the same (60 and 68 
seconds respectively). 

To summarize, 

• All independent events are reported as independent, with error prob- 
ability 8. 

• Ail first order pairwise event dependencies are found, with probability 
l-<?. 

• Some of the independent pairs reported may have dependency of higher 
order (undetectable by first order methods). 

As a negative example of the dependencies that can not be found by 
analysis of first order moments imagine that first event has a priori density 
of 2 events per maximal lag (AC) and is distributed uniformly everywhere 
except the IC secouds after occurrences of second event (which is, suppose, 
much sparsely distributed). Suppose also, that after each event of second 
type there are always two events of first type in time frame of K seconds and 
they are always separated by p or IC - p seconds. While the distributions 
of these two events are clearly dependent this fact can not be detected by 
analyzing the neighborhood event densities of first event around second event 
and vice versa. 



The presented methods can be extended to treat second-order depen- 
dencies, i.e., second-order moments of the distribution of one event in the 
neighborhood of another. One should be careful, however, to ensure that 
he has enough data to make a sufficient estimation of the measured quan- 
tities. This may be possible in the independent component analysis of the 
usual, non-sequential data, e.g., market basket data. In a market basket 
data all 'events' (purchases) in the dataset happen in the same time and 
we have multiple observations of the same variables (the set of all available 
goods). Removing the additional time variable may permit investigation of 
the higher order dependencies by approximating higher order moments witli 
sufficient precision. It also may allow analysis of more complex, non-pairwise 
dependencies. 



Table Sr. Experimental results with synthetic data. Found dependencies 
shown for each pair wise dependency detection method. 

Even t dependence analysis Uniform density Burst-like density 

Uniform event generation model (2, 3), (2,4) Every pair. 

Yes/No burst information (2,3), (2,4) Almost every pair. 

Burst density approximation (2,3), (2,4) (2,3), (2,4), (1,2). 

Episodes detection (3 -» 2), (4 -> 2) (3 ^ 2), (4 -> 2). 



In a preferred embodiment the underlying event -generation 
process is stationary, i.e., does not change over time, 
ergodic, and a quasi -Markovian property of the observed 
process; i.e., the distribution of events in the particular 
time frame depends only on some finite neighborhood of this 
frame. 

It will be appreciated that considerable advantages are 
achieved with the present invention when compared with 
solutions of prior art. A system and method operating 
according to features of the present invention, as 
contrasted to known methods, avoids a full search in the 
space of slices of event sequences by doing first a 
search in the k-dimensional Euclidian space. This 
approach has low computational complexity, which is 
linear in the size of the data. It will also be 
recognized, that the amount of memory storage during the 
prior art edit computations can be made smaller by using 
features of the invention. In addition, the computational 
time will be less using the method described in the 
present invention. Also, with lower computational time 
demands, simpler and cheaper processors can be utilized 
with the present invention but still achieving reasonable 
processing speed. 

It should be understood that the foregoing description is 
only illustrative of the invention. Various alternatives 
and modifications can be devised by those skilled in the 
art without departing from the invention. Accordingly, 
the present invention is intended to embrace all such 
alternatives, modifications and variances that fall 
within the scope of the appended claims. 
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